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I.  TJjJRODUCTION 

In  studies  of  free  electron  lasers  (FEL)  systems,  in  which  emission  gain 
is  induced  through  wiggler  magnetic  field  modulation  of  the  electron  motion, 
detailed  knowledge  of  the  electron  orbits  is  necessary  in  order  to  evaluate 
the  critical  features  of  the  amplification  scheme.  These  features  include 
evaluating  the  frequency  range  for  amplification,  gain,  rates,  saturation 
mechanisms,  and  amplifying  efficiencies.  However,  there  exist  very  few  mag¬ 
netic  geometries  where  the  electron  motion,  even  in  the  vacuum  fields,  can  be 
obtained  analytically  with  sufficient  information  to  assess  these  critical 


features.  Only  recently 


th  suffici 
l<i^has  it 


been  realized  that  in  the  "on  axis"  approxi¬ 


mation  the  helical  wiggler  orbit  problem  can  be  reduced  to  quadrature,  and  that 
particular  discovery  has  already  led  to  significant  progress  in  FEL  amplifica¬ 
tion  schemes. For  FEL  systems  with  strong  magnetic  guide  fields,  little 
or  no  Insightful  progress  has  been  made  except  for  this  helical  case.  However, 
for  just  this  class  of  magnetic  systems,  i.e.  strong  guide  fields  with  a 
relatively  low  amplitude  magnetic  wiggler  superimposed,  an  asymptotic  mathe¬ 
matical  formalism  does  exist  where  the  electron  equations  of  motion  in  the 
vacuum  fields  can  be  greatly  simplified,  and  typically  reduced  to  quadrature 
in  parameter  regimes  of  relevance  to  applications. 

For  wiggler  fields  which  are  small  in  amplf^ucte  compared  to  the  basic 
and  necessarily  geometrically  simple  guide  field,  asymptotic  methods  can  be 
used  to  simplify  the  orbits  for  various  ranges  of  the  parameter,  6=  kQv„/Qc  , 
where  fl  is  the  relativistic  electron  gyrofrequency  in  the  guide  field,  and 
kQ  and  v„  are  the  values  of  the  wiggler  wavenumber  and  electron  velocity 
parallel  to  this  guide  field.  From  the  current  studies  of  the  exactly  soluble 
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"on  axis"  helical  wiggler  these  same  ranges  are  in  fact  just  the  interesting  ones. 
For  6  «  1  ,  the  wiggler  madifications  of  the  quide  field  orbits  can  be  handled 
in  the  well  known  guiding  center  approximation,^  and  for  6  »1  ,  standard  weak 
perturbation  methods  easily  prove  to  be  adequate.  For  |1  -  5|  <  1  ,  the  so-called 
resonant  regime,  secular  perturbation  methods  have  proven  successful  in  simpli¬ 
fying  the  orbit  analysis  for  particle  motion  in  small  amplitude  fluctuating  rf 
fields.  This  analysis  can  be  straightforwardly  applied  to  the  wiggler  problem 
as  long  as  only  widely  separated  values  of  kQ  are  present  in  the  system,  and  as 
long  as  these  kQ's  are  very  slowly  varying  in  space.  Furthermore,  this  same 
method  can  be  used  when  6  is  close  to  the  ratio  of  two  discrete  integers.  How¬ 
ever,  of  particular  interest  to  FEL  applications  and  for  the  sake  of  simplicity, 
we  will  not  here  develop  a  general  formalism  for  an  arbitrary  wiggler,  but  will 
restrict  this  study  to  the  orbits  of  particles  near  the  magnetic  axis  of  a 
longitudinal  wiggler  field  in  cylindrical  geometry: 


1,2,3 


B  =  BqJ\  +  e  coskQz(l  + 


1  +  k?r2 


■4 


ez  +  2eokor  sink0zer 


The  exact  description  of  the  particle  motion  in  this  multiple  mirror  geometry 
can  be  greatly  simplified  because  of  the  existence  of  the  two  exact  constants 
of  the  motion--energy  and  angular  momentum--but  the  orbits  cannot  be  reduced 
to  quadrature  except  in  the  limit  6  -*■  0  Therefore,  the  success  of  the 
asymptotic  technique  for  this  case  is  particularly  meaningful.  Additionally, 
a  nunrt>er  of  FEL  studies  in  this  geometry  have  proceeded  neglecting  entirely 
the  radial  magnetic  field  in  Eq.  (1),  and  therefore  have  used  a  B  field  which 
is  not  divergence-free.  The  consequence  of  this  approximation  can  now  be 
determined. 
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II.  ASYMPTOTIC  TRAJECTORY  ANALYSIS 

For  the  specific  case  of  motion  In  a  longitudinal  wiggler  the  relativistic 
equations  of  motion  in  the  laboratory  frame  in  standard  cartesian  components 
are 

M  *  V  ■  -S-t,SV  -  W,<*Bx  *  14By>]  <2a> 

W  ■  S^vxSBy  '  VBx]  •  <2h> 

where  e  <  0  is  the  electron  change,  m  is  the  relativistic  mass  (which  is  a 

constant),  SB^  are  the  ith  vector  components  of  the  wiggler  field,  =  eBo/mc  , 

and  v  =  v  +  iv  .  The  entire  philosophy  of  the  trajectory  analysis  in  all  of 

x  y 

the  previously  discussed  ranges  of  6  is  that  the  left-hand  sides  of  Eq.  (2) 

induce  only  small  changes  in  the  orbits  described  by  the  right  side  of  Eq.  (2) 

in  a  cyclotron  period  T  *  2ir/nQ  .  With  such  a  philosophy  and  bv  Fourier 

analyzing  the  wiggler  fields,  a  straightforward  analysis  of  the  motion  can 

proceed,  even  in  the  resonance  regime.  However,  for  unification  and  simplicity 

o 

here  with  field  given  by  (1),  Eqs.  (2)  simply  become  to  order  e 


dv 

rft 


iftv  + 


1  to 


and 


(3a) 


dv  e«.krt 

sinkQz  tvuj*  -  v*u)j  ,  (3b) 

with  w  *  x  +  iy  ,  and  fl  =  nQ[l  +  c  coskQz(l  +  k^r2/4)] . 

Eqs.  (3)  are  not  exactly  Integrable,  so  a  perturbation  theory  must  be 
constructed  for  their  solution,  however  two  exact  constants  of  the  motion, 
energy  and  angular  momentum  (In  a  symmetric  multiple  mirror  system)  do  exist 
and  from  (3)  these  correspond  respectively  to  the  conditions 
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^<|v|z  *  v|)  -  0  (4a) 

obtained  by  multiplying  (3a)  by  v*,  and  adding  this  equation  to  its  complex 
conjugate  equation  and  using  (3b) 

and 


;  “*“  ♦  «m2]  =  o  (4b) 

[obtained  by  multiplying  (3a)  by  and  subtracting  this  equation  from  its 
complex  conjugate]. 

In  the  spirit  of  usual  secular  perturbation  theory,  the  solutions  of  (2) 
proceed  by  assuming  that  the  SB  corrections  induce  a  slowly  varying  component 
of  v,  and  a  slow  modulation  of  the  phase  and  amplitude  of  the  fast  component 
of  v.  So  writing 
v  =  v  +  6v 

(v,6v  are  the  slow  and  rapidly  varying  components  of  v  respectively)  and 
correspondingly 

(*)  *  u)  +  6u) 

v  =  V  +  6v  , 

2  2  Z 

the  slow  variation  of  v  is  determined  by 

£  *  V  •-£[«,*  -  v2(6Bx  ♦  ISBy)]  (5a) 

and  dv,  given  by 

(g^+  ifi0)6v  “  v2(5Bx  +  1<SBy^'  t1me  average)  (5b) 

and  similarly  for  vz,  and  Sv2,  where  the  bar  (  )  indicates  averaging  over  a 

cyclotron  period,  2ir/&0  .  Now,  the  quantity  S  *  k0vz/no  enters  in  directly 
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in  calculating  the  time  averages.  When  <5  «  1  (called  the  drift  approximate  or 
gyrokinetic  limit)  then  SB.  is  essentially  constant  in  a  cyclotron  period,  so 
here  we  have 


6  «  1 


3F  +  'V  me -  5z(5Bx  +  )] 


(6a) 


hence 


j  ev, 

5  ■  «  -  +  15V  • 


+  Sea) 2 


6B  i 6B  k  v 

*  y-r  *  V>  -  U'lriV  cosk„z(1  *  JS— L>  • 


(6b) 


We  see  from  (6b)  that  in  this  limit  v  consists  of  the  "x"  and  "y"  components 

of  the  velocity  parallel  to  B,  plus  the  usual  curvature  drift  term,  the  last 

term  in  (6b).  The  usual  grad  B  drift  is  neglected  in  this  PEL  application 
2  2 

because  |v|  «  v„,  so  it  is  correspondingly  smaller  than  the  curvature  drift. 
Even  so,  in  this  drift  limit,  the  curvature  drift  is  of  order  cfiku>0  times 
vz  ,  so  by  the  assumed  orderings  of  this  limit,  this  is  negligibly  small. 

Note  also  that  the  phase  of  the  curvature  drift  differs  from  id  by  a  factor 
of  ±i ,  hence  this  drift  is  in  the  "0"  direction,  i.e.  perpendicular  to  the 
vector  f  . 

Further  discussions  of  the  drift  limit  are  not  necessary  because  of  the 
vast  literature  on  the  subject.  However,  it  is  to  be  noted  that  an  overall 
amplitude  correction  is  to  be  added  (in  the  6  «  1  limit)  to  previous  orbit 
analyses7  which  entirely  neglected  6Bp .  This  correction,  easily  derivable 
from  Eq.  (2)  or  (3)  corresponds  to  use  of  the  adiabatic  invariant  p  *  e|6v|  /2B, 
and  corresponds  to  an  amplitude  factor  multiplying  the  previously  obtained 
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solutions^  for  u>  by  the  factor  This  geometric  factor  then  properly  takes 

into  account  that  the  actual  SB  must  satisfy  V*6B  =  0,  whereas  neglect  of  6Br 
does  not  satisfy  this  divergence  condition. 

1  2 

In  the  more  interesting  limit  for  obtaining  large  orbit  excursions,  * 

| 6 |  ~  1,  and  the  methods  of  time  averaging  the  exact  equations  must  be  care¬ 
fully  reconsidered.  In  particular  the  SB^  terms  vary  as  sinkQz  ,  tha  argument 
of  which  varies  in  time  roughly  as  kQvzt  +  slow  modulation.  Similarly,  from 
(5b)  we  see  that  6v  varies  as  exp(-iJ2Qt)*  a  function  with  a  slow  time  modulation. 
Therefore  in  time  averaging  (5a)  we  now  obtain,  to  lowest  order  in  e. 


i  efl  2  -  n 

u  «  -  - 2^-  exp(1*)[(l  ♦  +-M :>*  +  oo5/2i]  ,  (7b) 

and 

dv,  eJl  k  ,  . 

= - |_°.[(w  +  iv/no)  exp(it|»)u*  +  C.C.]  (7c) 


with  the  functions  6v  =  u  exp(-ifi  t),  -  k  z  +  Rt,  and  the  approximations 

0  0  0 

6u>  *  6v/(-in0)  [1  +  Ot\  \  i|>  «  Qq  ,  6vz  «  |6v|  (by  energy  conservation  6vz  ~  v6v/vz), 
and  it  was  assumed  that  k0w  remains  small  compared  to  unity.  As  v  varies  slowly 
compared  to  nQ,  (7a)  can  be  iteratively  solved  to  yield  (correct  to  order  e) 

.  (a,) 


u  - 


^■■T*«rnr  v.n 


Defining  p  =  exp(-1e2fl0t/8)u  and  e  «  -  e^/St  we  find  the  simple  equation  set 


eu. 

p  =  -  exp(i£)w 

v  *  £  *-|pexp(-1c)|l  ♦  fk+  ^  j 


d.r5  .  lo. 
dtlvz  2R0 


IpI2]  - 


Equation  (9c)  is  also  a  direct  result  of  kinetic  energy  conservation,  Eq.  (4a). 

Using  (9a)  in  the  lefthand  side  of  (9b),  and  neglecting  the  product  of  small 

2  2  2 

terms  times  e  (which  then  eliminates  the  kQr  /4  term  in  Bz,  but  is  valid  if 
kQr  remains  small) 

2  e2Q2 

P  -  i[k0vz  +  n0(i  -  ‘  1F  P  =  0  ‘  {10) 

This  equation  (10)  appears  to  be  as  complicated  as  the  original  starting 
equations  in  (3),  but  in  fact  now  there  is  only  one  time  scale;  the  slow  time 
scale  compared  to  2ir/fl0,  and  with  the  conservation  equation  (9c),  (10)  will 
be  shown  to  be  exactly  soluble. 

Introducing  the  new  dependent  variables  s  and  0,  p  *  flQs  exp(-fe), 

Eqs.  (8)  reduce  to 

2  q  s2  ft  k2  - 

8S2  ♦  .  c2  ,  (11a) 


M 

Tf 


£ 


C2  a  constant,  with  qQ  ■  k0Cj  +  (1  -  eV8)fl0  ,  Cj  the  constant  defined  by 
Eq.  (9c),  and 


s‘*“3sv(s>  * 


•it?  .he  effective  potential  V  given  by 


(lib) 


v  V<1o- 

2$*  0 


9  A  2  4 

2  nftk„qfts  fLK.  c 

x  o  0  0  X  o  o.o 


(lie) 
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2 

Choosing  the  dimensionless  variables  X  =  kQs  ,  n  =  qo/ftQ  ♦  and  ?  =  ^o^o  * 
we  have  the  final  equation,  soluble  by  quadratures 


By  their  definition  and  the  assumptions  of  this  analysis 

2  2  4  4 

X  =  ks«l,  n~l  +  <$«l,  and  t,  ~  s  k  n  +  k  s  «  1  , 

0  0 

therefore  the  time  scale  for  the  X  variation  is  seen  to  be  much  longer  than  Q0,  con- 

4  6 

sistent  with  assumptions.  We  note  that  the  X  and  X  terms  in  the  effective  potential, 
V,  come  strictly  from  the  change  in  vz  as  |z|2  increases  as  given  by  (9c)  and  is  a 
result  of  the  detuning  of  the  linear  resonance,  kQvz  «  .  Also  note  that  as 

<  0  for  electrons,  the  conservation  law  indicates  that  as  |z|  ( 6 v (  increases 

vz  decreases,  as  expected  from  constancy  of  y  =  [1  -  v2/c2]_Js.  The  X6  term  in 
the  potential  is  definitely  a  stabilizing  term,  that  is,  it  sets  an  upper  bound 
on  the  limits  of  |X|,  and  these  limits  are  not  such  as  to  invalidate  the  assump¬ 
tion  of  small  X  because  it  clearly  dominates  the  force  term  when  X  5  n  and  as 
n  «  1,  the  small  X  approximation  must  still  retain  the  X^  term. 

•  o 

By  conservation  of  the  effective  energy,  E,  of  Eq.  (12):  E  =  X  /2  +  V 

=  constant,  Eq.  (12)  can  be  reduced  to  quadratures  and  by  using  the  independent 
2 

variable  y  *  X  ,  the  integrals  are  all  of  the  elliptic  type.  However,  this 
result  is  probably  no  more  useful  than  the  observation  that  by  standardizing 
the  original  equation,  (3a),  and  neglecting  the  kQr  /4  term  in  6Bz,  the  resulting 
equation  is  of  the  Mathieu  form.  So,  some  general  considerations  are  in  order 
and  will  be  shown  to  provide  sufficient  information  for  analytically  describing 
the  electron  orbits  to  the  degree  necessary  for  most  FEL  applications. 
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III.  RESONANT  ORBIT  DESCRIPTION 

The  effective  potential  described  by  (12b)  yields  a  very  useful  method 
for  analytically  dealing  with  the  orbits.  Recalling  the  definitions  of 
variables  we  have  that. 

kQ6u)  =  IX  exp{-i  [Uo(l-e2/8)t+0]}  ,  (13a) 


so  that  the  effective  potential  diagram  describes  the  modulated  Lamor  motion 

of  the  electron  with  X  =  |k06oj|.  Knowing  the  quantitative  features  we  can 

derive  all  of  the  remaining  orbit  characteristics  including  the  guiding  center 
position, 

or 


„  r  <  e  n^n  o 

-~2-|e  -  -jjV  +  — 3?J=  constant  +  x 


(13c) 


Equation  (13c)  is  also  a  consequence  of  angular  momentum  conservation,  equivalent 
to  Eq.  (4b).  Therefore  following  the  X  motion  using  the  effective  potential  method 
directly  gives  a  picture  of  the  excursions  in  Larmor  radius,  |ko6u>(t)|  =  X,  and  the 
difference  in  X2  gives  the  differences  in  the  square  of  the  guiding  center  position; 

|ko<L(t1)|2  -  |kQ(L(t2)|2  =  X2(tl)  -  X2(t2)  . 

In  general  the  effective  potential  has  the  characteristic  shape  shown  in 
Fig.  1.  The  number  of  minima  can  either  be  one  or  two.  These  differences  are 
important  as  a  minimum  in  the  potential  corresponds  to  a  stable  electron  orbit 
with  constant  radial  guiding  center  and  magnitude  of  the  Larmor  radius.  The 
details  of  this  potential  and  the  resulting  motion  can  be  investigated  in  depth 
in  various  limiting  cases  of  the  parameter,  c,  e,  and  n. 

In  the  experimentally  interesting  limit  where  |n|  »  e/2,  the  equations  of 
motion  can  be  written  in  terms  of  the  new  variable,  y  =  X/ Ini'5,  and  new  potential. 


V.  as 


3? 

y  =- 3y 


y 


(14) 


J 


with  £(n)  =  {*j*  p  <  o'  c  ~  r>(kr)2n,  the  constant  c/n^  ~  (kr)2/n  and  is 
typically  of  order  unity.  From  (14)  we  see  that  the  time  scale  of  these  defec¬ 
tions  are  of  the  order  =2-rr/nQ0  »  2tt/S20  by  assumption.  By  simple  algebraic 

calculations  the  potential  V  can  be  shown  to  have  two  minima  (as  in  Fig.  1) 

2  2 

only  if  -2  <  <;/n  <  0,  or  n  <0  and  0  i  c/n  <  2/3  .  Otherwise  V  has  only 

one  minimum  for  y  >  0. 

As  a  particularly  simple  example  let  us  determine  the  quantitative  features 
of  a  particle  with  r,  =  0,  jn|  »  c/2.  For  n  >  0,  only  y  =  0  is  a  minimum,  and 
the  motion  is  oscillatory  but  y  remains  near  y  =  0.  For  n  <  0  two  minima  which 
are  zeros  exist  in  V  ,  at  y  =  0,  and  y  =  /8.  A  relative  maximum  of  v  occurs 
at  ym  =  /8/3  with  a  value  of  V(y^(  =  (16/27)n2ft2.  As  this  value  is  large 
compared  to  the  energy,  given  by  (13c)  as 

e2ft2  ?  2  2 

E  =  -^[Inwi  -  X2]  <  e2  , 

by  assumption  then  electrons  initially  near  the  minimum  at  yQ  =  /8  will  oscillate 

but  remain  near  to  this  value  of  y  .  For  small  oscillations  near  y0  the  motion 

is  sinusoidal  with  the  exact  frequency  being  given  by  uQ$  =  /+$  =  /Zhn^0  in  this 

"not  too"  resonant  regime,  |n|  »e/2,  one  or  two  equilibrium  positions  can  exist, 

and  motion  around  these  stable  singular  points  is  very  limited.  Additionally, 

in  this  not  too  resonant  regime  the  new  non-zero  guiding  center  position  is  of 

the  order  of  the  size  of  the  Larmor  radius,  and  is  n-dependent  (eg.  c  =  0  gives 

the  second  equilibrium  positions  as  |kQp|  *  /8|n| ) •  Hence,  small  changes  in  the 

resonance,  Ik  V  /Si  |  *  |6|  can  make  substantial  relative  changes  in  the  equilib- 
1  o  zo  o 

rium  orbit  positions. 


-12- 


2  2 

For  smaller  values  of  n>  very  near  resonance  such  that  n  -  e  /4  -  5  is 

negative  the  first  minimum  in  V  now  is  not  at  X  =  0  even  when  z,  -  0.  This  means 

that  even  those  electrons  that  initially  have  small  perpendicular  velocities  and 

radial  positions  near  r  =  0  can  still  experience  significant  radial  excursions 

and  velocity  modulations.  For  example  in  the  simplest  limit  c  =  0,  V  has  an 

absolute  minimum  at  X  *  X  .  with 

0 

X2/4  =  -5n/9  +(4/9)/nZ  +  (9/64)e2  (15) 

Therefore  an  electron  will  make  excursions  in  X  around  XQ,  turning  wherever  E  =  V(X). 
From  (13c),  we  note  that  for  £  =  0,  the  energy  E  =  (eZQZ/32)[X2  -  |kow|2].  So  for  a 
particle  with  an  initial  Larmor  radius,  (X),  equal  to  its  initial  guiding  center 

p  0 

displacement,  ( | kQu> | ) ,  E  =  0.  For  these  electrons  (c  =  0,  E  =  0,  0  <  e  /4)  the 

turning  points,  X^.,  can  be  shown  to  be  Xt  -  0  and  Xt  =  2(e/2-n)2.  Of  course  for 

the  particular  initial  condition  of  the  particle  initially  being  at  r  =  0  with 

exactly  zero  perpendicular  velocity,  i.e.,  add  the  special  condition  X  =  0  to 

these  initial  conditions;  the  particle  just  remains  at  X  =  0.  However,  this  is 

an  unstable  point  and  very  small  deviations  from  the  conditions  X  =  |k0w|  =  0 

will  cause  the  electron  to  make  significant  radial  excursions. 

In  the  simplest  case  of  z,  -  0  we  should  note  that  the  potential  has  two 
2  2 

minima  if  n  >  e  /4  and  q  <  0  and  these  minima  are  at  X  =  0  and  X  =  XQ.  If 
E  =  0  in  this  circumstance  the  particles  stability  remain  at  X  =  0  and  X  =  XQ. 

This  demonstrates  another  parameter  range  where  electrons  injected  with  only 
a  small  perpendicular  velocity,  but  on  the  r  *  0  beam  axis  will  not  remain  there. 
In  this  case  of  E  =  0,  the  electron  would  sweep  through  r  =  0  every  Larmor  cycle 
with  Its  guiding  center  remaining  at  the  fixed  distance  XQ  =  |k0w|*  from  the  axis, 
and  slowly  rotating  in  angle. 
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In  both  of  these  cases  with  n  -  e  /4  negative  or  positive,  XQ  or  the 
turning  point,  X*,  defines  the  amount  of  excursion  possible:  Xmav  »  X.  =  X 
-  j kQuj | ,  -|k05vA/^|,  and  as  such  measures  both  the  perpendicular  "thermal" 
spread  that  a  beam  can  develop;  and  its  typical  guiding  center  displacement. 
For  example  when  n  «  e,  Xt  *  /2 c  ,  giving  a  maximum  radial  displacement  of  an 
electron  of  order  2/2e/kQ;  and  a  perpendicular  velocity  6vA  =/2evz.  Both  of 
these  displacements  are  significant  and  must  be  treated  carefully  and  self- 
consistently  in  the  modeling  of  electron  beam  equilibria  and  free  electron 
gain  factors. 
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IV.  DISCUSSION 

The  analysis  presented  here  demonstrates  that  in  the  vicinity  of  resonance, 

jkzvz/fi0|  *  1  analytic  techniques  can  be  used  to  obtain  the  electron  orbits 

even  in  a  longitudinal  wiggler  including  the  radial  component  of  the  magnetic 

field.  It  is  also  possible  to  allow  for  adiabatic  variation  of  the  guide  field 

and  wiggler  wavelength,  within  this  formalism  if  the  quantities  vary  slowly  in 

a  period  of  this  reduced  motion,  which  is  of  order  2n /Q  e,  or  2ir/nft  whichever 

o  o 

is  shortest.  This  same  type  of  secular  analysis  has  been  applied  to  both  the 
helical  and  linear  wiggler  geometries  without  the  assumption  used  here,  of  the 

g 

electrons  being  near  the  magnetic  axis  of  the  guide  field.  Additionally,  a  more 
general  perturbation  formalism  has  been  applied  to  a  general  wiggler  geometry 
with  simple  periodicity  in  Z  assuming  only  a  small  wiggler  field  strength  to 
guide  field  strength.  The  question  of  stochasticity  of  the  electron  motion  has 

O 

also  been  addressed. 

In  all,  very  important  questions  of  the  quality  of  the  in  situ  electron 
beam  can  be  answered.  Equally  significant  is  that  in  the  vicinity  of  resonance 
|k»v„/ft0|  =  N,  N  integer,  where  large  radiation  conversion  efficiencies  can  be 
expected,  analytic  approximations  to  the  electron  orbits  can  be  obtained  and 
these  are  all  important  in  determining  realistic  gain  factors  and  saturation 
mechanisms  and  radiation  strengths  for  actual  free  electron  lasers. 
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